The local content of bipartite qubit correlations 
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One of the last open problems concerning two qubits in a pure state is to find the exact local content of their 
correlation, in the sense of Elitzur, Popescu and Rohrlich (EPR2) [Phys. Lett. A 162, 25 (1992)]. We propose 
a new EPR2 decomposition that allows us to prove, for a wide range of states \ip(9)) = cos 9\00) + sin 
that their local content is p~L (6) = cos 28. We also share reflections on how to possibly extend our result to all 
two-qubit pure states. 



INTRODUCTION 



The incompatibility of quantum mechanics with local vari- 
able theories, as shown by Bell [1], lies at the statistical level: 
local variable theories cannot reproduce all statistical predic- 
tions of quantum theory. In a typical Bell experiment, one 
observes correlations between the measurement results of two 
partners (Alice and Bob), and averages them over measure- 
ments on many pairs of particles. One may then conclude that 
non-locality was observed if the average correlations thus ob- 
tained violate a Bell inequality. 

However, if the statistics of the observations exhibit non- 
locality, it does not imply that all individual pairs behave 
non-locally. This observation lead Elitzur, Popescu and 
Rohrlich |2|], hereafter referred to as EPR2, to wonder whether 
one could consider that a fraction of the pairs still behaves lo- 
cally, while another fraction would behave non-locally (and 
possibly more non-locally than quantum mechanics allows). 

More explicitly, writing Pq the quantum mechanical prob- 
ability distribution for Alice and Bob's results, the EPR2 ap- 
proach consists in decomposing Pq as a convex sum of a local 
part, Pl, and of a non-local part, Pjvx, in the form 



THE EPR2 APPROACH FOR TWO-QUBIT PURE STATES 
Correlations of two-qubit pure states 

Without loss of generality, any two-qubit pure state can be 
written in the form 



\tp{6)) = cos(9|00) + sin0|ll) 



(2) 



with 9 G [0, ?]. In the following, we shall use the notation 
c = cos 29, s = sin 29 (with c, s G [0, 1]). 

Each qubit is subjected to a Von Neumann measurement, 
labeled by unit vectors a and b on the Bloch sphere S 2 . Let 
us denote by a z and b z the z components of a and b, by 
aj_ = \Jl — a 2 and b^_ = J\ — b\ the amplitudes of the 
components of a and b in the xy plane, and by x G ] — 7r, 7r] 
the difference between the azimuthal angles of a and b 1 , with 
b 1 defined as the reflection of b with respect to the xz plane 1 . 

With these notations, and for binary results a, (3 = ±1, 
quantum mechanics predicts the following conditional proba- 
bility distribution: 



P Q (a,f3\a,b) 



aM Q (a) + (3M Q (b) 
+ a/3E Q (a,b)) 



(3) 



X 
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Pq=PlPl + (1-Pl)Pnl, With p L G [0,1]. (1) 



The maximal weight p~L = maxpi that can be attributed to 
the local part can be regarded as a measure of (non-)locality of 
the quantum distribution Pq. Finding this maximal possible 
local weight is not a trivial problem; so far one only knows 
how to calculate lower and upper bounds on p~L. 

In this paper, we concentrate on the simplest case of a quan- 
tum probability distribution originating from Von Neumann 
measurements on two-qubit pure states. After recalling pre- 
viously known results for this case, we propose a new EPR2 
decomposition and derive a new lower bound on p~L, which 
reaches the previously known upper bound Q| for a wide class 
of states. This gives a definite value for the exact local con- 
tent p~L of those states. We then share reflections on how one 
might possibly extend our result to all two-qubit pure states. 



with Mq(o) = ca z , M Q (b) = cb z , 
E Q (a, b) = a z b z + s (a x b x - a y b y ) 



a z b z 



cos x ■ 



(4) 
(5) 



As explained before, the EPR2 problem is to find a decom- 
position of -Pq(«, f3\ a, b) as a convex sum of a local and a 
non-local probability distribution, in the form (HJ. For a given 
state (i.e., a given value of 9), the equality is required to hold 
for all possible measurements a, b and for all results a, (3. The 
weight pl G [0, 1] of the local distribution should be indepen- 
dent of the measurements and of the outcomes. 



1 The axes x,y,z of the Bloch sphere are defined as usual: |0) and |1) 
are identified with the north and south poles (i.e., along the z axis), 
while the state (|0) + |l))/v2 defines the x direction. We introduce 
b' = (b x , —by, b z ) to account for the minus sign in front of a y b y in (f5j- 
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The probability distribution Pl(&, (3\ a, b) is required to be 
local, in the sense that it can be explained by local variables 
A, i.e. it can be decomposed in the form 



P L (a,P\3, b) 



d\p(\)Pf(a\a)pZ(p\b). 



(6) 



On the other hand, no restriction is imposed on Pnl, except 
that it must be non-negative for all inputs and outputs: 



Pnl = 



1 



l-VL 



(Pq-PlPl) > 



(7) 



In particular, Pn l is allowed to be even more non-local than 
quantum mechanical correlations 2 . 

The goal is to find, for a given state, a decomposition with 
the largest possible value for pl, denoted p~l(9), which char- 
acterizes the locality of the probability distribution Pq as de- 
fined by OHD. 



Previously known results and conjecture 

In their original paper J2l, Elitzur, Popescu and Rohrlich 
proposed an explicit local probability distribution Pl, which 
lead to an EPR2 decomposition with pl = t- 2 - This was the 
first known lower bound on P~l(6). Clearly, this was not opti- 
mal, at least when approaching the product state (9 = 0, i.e., 
s = 0) which is fully local, and therefore satisfies p~l(0) = 1. 

They also argued that for the maximally entangled state 
(9 = 7r/4), Pq contains no local part: p~l(j) = 0, i.e. no 
EPR2 decomposition with Pl > exists for this state. This is 
in fact a much more general result, as shown later by Barrett 
et al J3]: the maximally entangled state of two e?-dimensional 
quantum systems, for any dimension d, has no local compo- 
nent. 

In J3], one of the authors could improve on the first lower 
bound for p~l(0), as he gave an explicit decomposition that 
achieves pl = 1 — s. Interestingly, it was noted that this 
is the largest possible value that can be attributed to pl, if 
Pl depends only on the z components a z and b z of a and b. 
Recently, this bound (pi > 1 — s) has been extended to mixed 
two-qubit states by noticing that s is actually the concurrence 
of the state ]5|] . 

It is worth mentioning here that finding lower bounds on 
p~L (ie., explicit EPR2 decompositions) can be useful for the 
problem of simulating quantum correlations with non-local 
resources, as only the non-local part then needs to be simu- 
lated. The decomposition of ^ was thus successfully used to 
simulate partially entangled 2-qubit states @]. 



On the other hand, an upper bound on p~l(0) can be ob- 
tained with the help of Bell inequalities J3]. Let / < II be a 
Bell inequality (defined by a linear combination of conditional 
probabilities), Iq the quantum value obtainable with the prob- 
ability distribution Pq, and Ins(> II) the maximum value 
obtainable with non-signaling distributions. Then from (Q~|i it 
follows that Iq < p L I L + (1 - pl)Ins, i.e., 



PL 



< 



Ins ~ Iq 



(8) 



LNS 



Using the family of "chained Bell inequalities" QVL|8|] , an up- 
per bound for p~L (9) was derived (numerically) in yfl, namely 

p- L {9) < cos 29. 



So far, the gap was still open between the two bounds 

1 - sin 29 < p L (9) < cos 29 . 



(9) 



It has been conjectured |@] that there should exist an EPR2 
decomposition that reaches the upper bound, i.e. withp^ = c. 
If this could be proven to be true, then the lower and upper 
bounds would coincide, and one could conclude that the value 
of p~l(9) is exactly cos 29. 

In the following we describe our (partially successful) at- 
tempts to prove this conjecture. 



REFORMULATION OF THE PROBLEM 
TO PROVE THE CONJECTURE 

Our goal is now to see whether it is indeed possible to at- 
tribute a weight pi = c in the EPR2 decomposition of the 
2-qubit probability distribution Pq (0, and write 



Pq = cP L + (1 - c)P, 



NL 



(10) 



For that, we want to find an explicit local probability distribu- 
tion Pl, such that Pnl = jzt^.{Pq — cPl) is a valid probabil- 
ity distribution, i.e. Pnl must be non-negative. The problem 
thus translates into 

Problem : find P L , such that P Q - cP L > . (11) 



At this point, we shall impose an additional (and possi- 
bly questionable) constraint on the EPR2 decomposition we 
are looking for: we want the non-local part to have ran- 
dom marginals 3 , i.e., with obvious notations, Mjvi(a) = 
MnlQ}) = 0. The intuition is that the marginals are local 



2 Note however, that Pjvi = i— PL (Pq ~ PlPl) is by construction non- 
signaling. 



3 Note that this constraint precisely justifies the choice pl = c. Indeed, 
if one can find an EPR2 decomposition with random non-local marginals, 
then for the setting z, Mq(z) = c = plMl{z), which implies pl > c. 
Now, c is known to be an upper bound for pl , and therefore Pl = c- 
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properties, which should be concentrated on the local compo- 
nent only 4 . 

As equality ( fTOb should also hold individually for the 
marginals on Alice's and Bob's sides, one should then have 
M Q (a) = cM L (a) and M Q (b) = cM L (b), i.e. 

M L (a)=a z , M L {b)=b z . (12) 
With these constraints, the condition Pq — cPl > reads: 

for all a, f3, a, b, 

l-c + aP(E Q (a,b)-cE L {a,b))>0. (13) 

Thus, the problem now translates into: 

Problem : find Pl , such that 
M L (a)=a z 

M L (b)=b z (14) 
\E Q (a,b) - cE L (a,b)\ < 1-c 

PROPOSAL FOR A NEW EPR2 DECOMPOSITION 

As we are dealing with qubits, the natural geometry of the 
problem involves unit vectors on the Bloch sphere; we shall 
propose a local component Pj, that makes the most of this 
geometry. Inspired also by models that Bell devised to repro- 
duce the measurement statistics on a single qubit in the state 
|0) Jl2ll (which gives precisely the marginals we want), or to 
approximate the statistics of the singlet state |l|], we introduce 
the following model to define Pl : 

Local model: Alice and Bob share a random local vari- 
able A, uniformly distributed on the Bloch sphere. When Alice 
receives the measurement direction a, she outputs ct(a, A) = 
sign(a z — a ■ A). Similarly, when Bob receives the measure- 
ment direction b, he outputs [3{b 1 A) = &ign{b z — b' ■ A), where 
b' is the reflection ofb with respect to the xz plane. 

Let us check whether the constraints ( fT4T > are satisfied. 

Marginals. Alice's and Bob's marginals corresponding to 
our local probability distribution Pl are, as required in (14\ : 



M L (a) 



dA 

4.7T 

6X 



sign(a z - a • A) = a z (15) 



M L (b) = Ij — sign(b z - b' ■ A) = b z . (16) 



4 However, one hint that the argument is questionable is the fact that the 
no-signaling polytopes for arbitrarily many measurements but binary out- 
comes contain extremal points with non-random marginals HolfTlll . 



Correlation term. The details for the calculation of the lo- 
cal correlation coefficient El (a, b) are given in Appendix A. 
We find 



E L (a,b) = 



^ sign(a 2 — a - A) sign(6 z - b' • A) 

S 2 47T 



l-|a,-6,| ifx = 0, 
\a z + b z \-l ifx = 7r 



1 - 



2|X| 



' a±b z —a z b± cos 



(17) 



t 1 sO,arctan(^ b±sinlxl 
-H z arctan( a '^- a ^ fc r i° sx ) 

7r A ^ aj_ sin Ixl J 

if < Ixl < 7T . 



One can then check 5 that for all settings a and 6, 

\E Q (a, b) - cE L {a, 6)| < max(l - c, c - s) . (19) 

The last constraint in ( TPfl i is thus satisfied when c — s < 
1 — c, i.e. when c < |: 

For all c < 0.8, for all o, b, 

\E Q (a,b)-cE L (a,b)\<l-c. (20) 



CONCLUSION REGARDING OUR EPR2 
DECOMPOSITION 

When c < 0.8, since our local probability distribution P^ 
satisfies the three constraints (IT4t . it defines a valid EPR2 de- 
composition for Pq, with a local weight that can take the value 
Pl = c. This gives the lower bound p~L (6) > cos 28 for all 
pure two-qubit states (|2]) such that cos 29 < 0.8 (or 9 > O.Itt). 
As cos 29 is also known to be an upper bound for pl {9) 13J, 
we conclude that this is actually its definite value: 



when cos 29 < 0.8, p~ L (9) = cos 29 . 



(21) 



When cos 29 > 0.8 however, there exists measurement set- 
tings a, b for which the third constraint in (fT4b is not satisfied 
by Pl 6 - Our local probability distribution cannot be attributed 
a weight pl = c in that case. 

Still, our decomposition gives a non-trivial lower bound 
on pl even when c > 0.8, namely 7 pl > c + s — 1 + 



5 This can be proven analytically for the cases when \ = or it: for a 
given c, one gets the result by looking at the maximum of the function 
\Eq — cEl\ for all a z ,b z . For < |x| < T on the other hand, we 
checked the bound i 1 9t numerically; for each value of c, it was only 3 
parameters to vary, so we are confident that the numerics are trustworthy. 

6 Take a = —b = x for instance: Eq(x, —x) — cEl{x, —x) = c — s > 
1 - cif c> 0.8. 

7 The lower bound is given by min „ _. r -5^-. As for the case c < 0.8, 
the bound can be obtained analytically for \ = or 7r, and was checked 
numerically for < |x| < 7r, 
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Upper bound: p < c 




0.05 0.1 0.15 0.2 0.25 

6/JI 



FIG. 1 : Bounds on ph (6) in the EPR2 decomposition for two-qubit 
pure states, as a function of the state parameter 6. Our new lower 
bound pl (6) > cos 26, valid for all states such that cos 26 < 0.8, 
coincides with the upper bound obtained in (J. There is still a gap 
between the lower and upper bound when cos 29 > 0.8. 



extreme cases, we would still impose A a = Ab for the max- 
imally entangled state (9 = j), while the two A's would be 
completely decorrelated for the product state (9 = 0). 

The problem is now to find the proper way to correlate the 
two A's for each state, i.e. determine the distribution functions 
pe(X a , Ab). Here are a few properties that we might want to 
impose on pe(X a , X b ): 

• Forgetting about A;,, A a should be uniformly dis- 
tributed, and vice versa. This will ensure in particu- 
lar that the marginals are those expected: Mj,(o) = 
a z , Ml (6) = b z . One should thus have: 

for all A a , JJ dX b p 9 (X a , Xb) = ^- ; (22) 
for all A 6 , JJ^dX aPs (X a , X b ) = ^ . (23) 



y/2(l - c)(l - s). As long as c < § (or 9 > 0.06tt), 
this lower bound is larger than the previously known bound 
1 - s II, but when c > y| our new decomposition gives a 
smaller bound. 

Figure[T]summarizes all the bounds we now know on pl(9). 
PROSPECTS 

We thus could prove the conjecture that pi (9) = cos 29 for 
all states such that cos 29 < 0.8. This reinforces our opinion, 
that the result should indeed hold for all pure two-qubit states. 

Unfortunately, we could not find so far an EPR2 decompo- 
sition with pl = cos 29 for the very partially entangled states 
(such that cos 29 > 0.8). Let us however share a few reflec- 
tions on how one could possibly look for a suitable local com- 
ponent Pl, which would allow one to prove the conjecture in 
full generality. 



• Let us denote by ($ a (b)>¥>a(b)) the spherical coordinates 
of A o ( ) . It looks very natural to impose that pg (X a , Xb) 
should only depend on <& a , $ b and ip = ip b — ip a , that 
it should be symmetrical when exchanging A a and A , 
and that it should have an even dependence on <p: 

Pe{X a ,X b ) = pe(tf a ,$ b ,<p) (24) 
= p 9 {p h ,# a ,<p) (25) 
= po(0 a ,&b,-<p) . (26) 

• According to an argument presented in Appendix B, not 
all pairs (A a ,A&) should be allowed. More precisely, 
writing c a = cos^,s a = sin 4r,c = cos ^Sb = 
sin ^ and c v = cos (p, one should have 

Pe{X ai X b ) = if -— ^ — < a . (27) 

-t CqC^C^ 



We realize that our local distribution Pl above fails to sat- 
isfy the constraints ( fPfl i when the state under consideration 
becomes less and less entangled. In our local model, it might 
be that we correlated the two parties too strongly, by imposing 
that they share the same local variable A. 

One idea would be to provide the two parties with two 
local variables A a and Ab G S 2 , while still considering re- 
sponse functions of the form a(a, A a ) = sign(a z — a- A a ) and 
f3(b, Ab) = sign(&£ — b 1 ■ X b ). Instead of imposing A a = Ab 
as in our previous model, we would correlate A a and Ab in a 
smoother way, depending on the state we consider 8 . In the 



We therefore suggest the following research program, to 
prove the above conjecture for all states: find candidate func- 
tions pe(X a ,X b ) that have the previous desired properties (1231 - 
|27| |. and then check whether the induced local probability 
distributions Pl satisfy the constraints (fl4T >. If one can find 
such solutions, then this will prove that pl(6) = cos 29 in- 
deed holds for all two-qubit pure states. On the other hand, 
if it turned out to be impossible to find such a function, then 
we might need to change our local model, and maybe relax 
the assumption that the non-local part should have random 
marginals. 



To prove the conjecture for — * 0, it is actually necessary to have a local 
part that depends on the state, contrary to our first proposal. Indeed, in the 

first order in 8 (or s), the constraint \Eq — cEl\ < 1 — c implies that El = Eq + 0(8) = a z b z + sa±b± cos \ + o(s). 
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Appendix A: Calculation of El (a, 6) 

Here we calculate the correlation coefficient E L (a, b) for 
our local probability distribution P L : 



E L {a, b) 



g2 4ir 
dA 

g2 47T 



sign(a 2 — a ■ A) sign(6, — b' ■ A) 

(1 - 2[a ■ A > a z \) (1 - 2[b' • X > b z 



= a z +b z -l + - d\[a-\>a z ] [b'-X>b z ], 
7T JJs 2 

(28) 

where [•] is the logical value of what is inside the brakets. The 
integral represents the area of the intersection of two spherical 
caps centered around a and b', and tangent to the north pole 
of the Bloch sphere. 

Let us parameterize A e S 2 by its zenithal and azimuthal 
angles (■$, ip), where <p is defined (for simplicity) with respect 
to the vertical half-plane that contains a. As El (a, b) should 
not depend on the sign of x (the difference between the az- 
imuthal angles of a and b'), it is sufficient to calculate it for 
X > 0, and simply replace x by |x| m the final expression. 
Also, we assume for now that a and b are both in the north 
hemisphere of the sphere. 



The two spherical caps can then be defined as 

{A | a ■ A > a z } 

= {(^)l^eh|,|]andtfe[o,^„M]} 



{X\b'-X>b z } 
withCM>CM € [O) 71 "] such that 



= m <p) \<P e [X - \x + ~] and G [0, CM]} 



cos C(^) 



cos cm = 



2 2 2 

a z — a ± cos <p 
o? z + cos 2 (p ' 

b\ -b\ cos 2 (y-x) 
b 2 z + b\ cos 2 (<p — x) 



Let us define ipo as the azimuthal angle for which 
CM) = CM)- The integral in d28j) can then be cal- 
culated as follows: 



dA [a ■ A > a z ] [b' ■ X > b z ] 

7r/2 rniin(ti£.(<p),tiZM) 

dip / sin^di? 

X-t/2 ^0 

d^fl-cosCM] 

X-ir/2 
,tt/2 

+ / d^[l-cosCM] • 
J if 



2 2 
z~ a ± 

2 +a 2 ± cos^ ip 



2 2 2 

Using the antiderivative fdy ^f""- 1 - cos v 
2a z arctan(a z tan ip) — tp, we find: 

dA [a ■ A > a z ] [b' ■ A > b z \ 

= 2tt — 2\ + 2a z arctan(a z tan ipo) — ira z 
+2b z arctan(& 2 tan(x — <^o)) — n b z . 

(29) 

We note that CM) = ^m('Po) implies a z 6j_cos(x — 
ipo) = a±b z cos <po, which in turn implies a 2 taniy9o = 

a^b -a^ cos X an( j b t , - l p )= *:b ± ~a^b z cos X 

serting these values in ( |29l ), then inserting the integral in d28l ) 
and writing |x| instead of x, we get the correlation coefficient: 

tp f- t\ 1 2 IaI , 2 , a ± 6 z - a 2 6 ± cosx, 

EL(a,b) = 1 1 — a z arctanf ; ; — : ) 

7r 7T 0j_sm|x 

2 ,a z bj_ - a±b z cosx, 

H — b z arctanl : — : — : ) . 

7T aj_ sm |x| 

(30) 

So far we have calculated this coefficient for settings in the 
north hemisphere of the Bloch sphere. If the settings are in 
the south hemisphere, one can use the fact that E L (a, b) = 
—E L (a,—b) = —E L {—a,b) = El{— a, — b). One can check 
that the above expression is actually still valid for all cases. 
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Note finally that in the above calculation, we didn't pay at- 
tention to particular cases, when the denominators in the frac- 
tions would be zero. For x = or tt, El (as in eq ( fT8l )) can 
be obtained by taking the corresponding limit in the previous 
expression, or can be obtained directly in a much simpler way. 



Appendix B: Allowed pairs (A , \b) in our last proposal 

Here we will argue that in our last proposal with two differ- 
ent local variables A a and \ for Alice and Bob, not all pairs 
(A a , X b ) should be allowed. 

The argument is based on the following observa- 
tion: suppose that Alice measures along direction a = 
(aj_ cos</3 a , a± smifa, a z ), and finds the outcome a = +1; 

this projects Bob's state onto (+a\<S>l \ip{0)) = y / l+ ^ a ' \b s ), 

with b s = cos^ a , sin^ a , If Bob then 



measures the setting 6g, he will necessarily get the result 
[3 = +1, and therefore 

P Q (+-\a,b s )=0. 

This in turn implies, for the EPR2 decomposition Pq = 
PlPl + (1 - Pl)Pnl (with p L ^ 0), that 

P L {+- \a,b s ) = 0. 

This constraint must be satisfied by any setting a (which 
defines the setting bs). To ensure this, we shall not allow pairs 
(A a ,Ab) that may give the results (a = +l,/3 = —1), for 
some choice of settings of the form (a, bs). 

To make this more explicit, let us fix the first local variable 
A a . For simplicity, we assume that A a is in the xz plane. If 
this is not the case, the analysis below would be slightly more 
tedious, but the final result would be the same. 

The settings a that give the result a(a, X a ) = +1 span the 
half-sphere A above the bisector plane between z and A a ; see 
Figure |2] (left) for a 2D representation. A can be defined as 

A = {a\u ■ a > 0}, where u ~ (— c a , 0, s a ). 

(Let us recall the notations: (i? a (b);V 3 a(6)) are the spheri- 
cal coordinates of A q (m, and we write c a = cos%-,s a = 



cos ■ 



•-, s b = sin and c v = cos(tp b - ip a ).) 



The settings bs, corresponding to these settings a G A, then 
also span a spherical cap, £>, included in A 9 . Using the fact 



that Or = } bs \l ^ 
defined as 



and a± 



-, B can in turn be 



Ai,(A„, C ) 





FIG. 2: Construction of the set A(,(A a , c) of allowed local variables 
Ab, given the first variable A a . Left: 2D cut in the vertical plane that 
contains A a . Right: 3D representation of Aj (X a , c) . For both figures, 



d a = f,C=0.5. 



According to the above observation, the allowed local vari- 
ables Af, must be such that for all those settings bs in B, 
P(bs,X b ) ^ —1, i.e., (Ah — z) ■ b' s < 0. This implies that 
X b — z should be in a cone centered around v, and with a half- 
angle £ = arcsin This writes 

v ■ > | \v 1 1 cos f = s . 

\\Xb - z\\ 

For the fixed X a considered here, the set of allowed vari- 
ables X b is then the intersection of this cone, translated by 
z, and the Bloch sphere; see Figure [2] Writing Ab = 



(2c b .s b c v , 2c b s b s v , 
that: 



1 — 2sl), the previous condition implies, 



S a S b 



the pair (A a , Ab) should be allowed only if — > s . 

1 C a CbC^ 

This justifies the constraint ( l27l l. 



B = {bs\u -a>0} = {b\-v-b> cs a }, 
where v = (sc a ,0,~s a ). 



9 Note: in particular, if \ a is not assumed to be in the xz plane, one would 
have 6' C A. instead, where £>' is the reflection of B with respect to the 
xz plane. 



